Analysis of multipath interference in three-slit experiments* 
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It is demonstrated that the three-slit interference, as obtained from explicit solutions of Maxwell's equations 
for realistic models of three-slit devices, including an idealized version of the three-slit device used in a recent 
three-slit experiment with light (U. Sinha et aJ., Science 329, 418 (2010)), is nonzero. The hypothesis that the 
three-slit interference should be zero is the result of dropping the one-to-one correspondence between the sym- 
bols in the mathematical theory and the different experimental configurations, opening the route to conclusions 
that cannot be derived from the theory proper. It is also shown that under certain experimental conditions, this 
hypothesis is a good approximation. 
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, According to the working hypothesis (WH) of Refs. 11112], quantum interference between many different pathways is simply 
the sum of the effects from all pairs of pathways. In particular, application of the WH to a three-slit experiment yields yj 

/(r,000) = |v/i(r) + V/2(r) + r3(r)P, (1) 

where xifj with j = 1,2,3 represents the amplitude of the wave emanating from the jth slit with the other two slits closed and 
1^ r denotes the position in space. Here and in the following, we denote the intensity of light recorded in a three-slit experiment 
Ci by /(r,000), the triple O's indicating that all three slits are open. We write /(r,COO) for the intensity of light recorded in the 
experiment in which the first slit is closed, and so on. 
Assuming the WH to be correct, it follows that 
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I. INTRODUCTION 



/(r,000) = |v/i(r) + v/2(r) + r3(r)r 
^ ' = |v/i(r) + r2(r)P + |V/i(r) + r3(r)P + |v/2(r) + V/3(r)p-|v/i(r)p-|v/2(r)|2-|v/3(r)|2 

\0'. = /(r,OOC)+/(r,OCO)+/(r,COO)-/(r,OCC)-/(r,COC)-/(r,CCO). (2) 

' In other words, still assuming the WH to be correct, we must have 

A(r) = /(r,000)-/(r,OOC)-/(r,OCO)-/(r,COO)+/(r,OCC)+/(r,COC)+/(r,CCO) =0. (3) 

In analogy to the expression for the two-slit interference term /(r, OO) — /(r, OC) — /(r, CO) in a two-slit experiment we refer 
to A(r) as the three-slit interference term. 

According to Refs. [110, the identity Eq. (O follows from quantum theory and the assumption that the Born rule /(r) °^ |^(r) p 
holds. In a recent three-slit experiment with light |2], the seven contributions to A(r) were measured and taking into account the 
?H uncertainties intrinsic to these experiments, it was found that A(r) « 0. This finding was then taken as experimental evidence 
5r that the Born rule /(r) oc |^(r) p is not violated |2]. 

The purpose of the present paper is to draw attention to the fact that within Maxwell's theory or quantum theory, the premise 
that Eq. ([T]| (which implies Eq. (O) holds is false. By explicit solution of the Maxwell equations for several three-slit devices, 
including an idealized version of the three-slit device used in experiment fJ], we show that A(r) is nonzero. We also point 
out that summing up the seven contributions to A(r), being the outcomes of seven experiments with different slit configurations. 
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FIG. 1: (Color online) Amplitudes of the £, (left) and E, (right) components of the electric fields as obtained from a FDTD solution of 
Maxwell's equation for light incident on a metallic plate with three slits. The incident wave is monochromatic and has wavelength A = 500nm. 
The slits are X wide, their centres being separated by 3A. The index of refraction of the 4A -thick metallic plate (colored black) is 2.29 + 2.61i. 
In the FDTD simulations, the material (steel) is represented by a Drude model 0]. 




FIG. 2: (Color online) Left: Normalized angular distribution of light /(e;000)//(e = 0;OOO) transmitted by iV = 3 slits (see Fig.[T} as 
obtained from the FDTD simulation (bullets) and Fraunhofer theory (solid line) l{d,s = \,d = 2,N = 3), see Eq. l|6ll, where s and d are the 
dimensionless slit width and slit separation, respectively |4]. Right: Normalized difference = [/(0;OOO) — /(0;COO) — /(0;OCO) — 
/(e;OOC) +/(e;CCO) +/(e;COC) +/(e;OCC)]//(e = 0;OOO) as a function of 6. According to the WH of Refs. QUI], this difference 
should be zero. 



requires dropping the one-to-one correspondence between the symbols in the mathematical theory and the different experimental 
configurations, not satisfying one of the basic criteria of a proper mathematical description of a collection of experiments. How- 
ever, as we also show, under certain experimental conditions, Eq. (O might be a good approximation. We present a quantitative 
analysis of the approximative character of the WH Eq. ([T]i and discuss its limitations. 



II. SOLUTION OF MAXWELL'S EQUATION 



The approximative character of the WH Eq. ([T]i can be demonstrated by simply solving the Maxwell equations for a three-slit 
device in which slits can be opened or closed (simulation results for the device employed in the experiment reported in Ref|2 
are presented in SectionUVTl. For simplicity, we assume translational invariance in the direction along the long axis of the slits, 
effectively reducing the dimension of the computational problem by one. 



A. Computer simulation of a three-slit device 

In Fig. [T] we show the stationary-state solution of the Maxwell equations, as obtained from a finite-difference time-domain 
(FDTD) simulation |3] for a three-slit device, with slits being X wide and their centres being separated by 3A, illuminated by a 
monochromatic wave with wavelength X. From the simulation data, we extract the angular distribution /((J,000). Repeating 
these simulations with one and two of the slits closed, we obtain 7(0, COO) and so on. In all these simulations, the number of 
mesh points per wavelength A was taken to be 100 to ensure that the discretization eiTors of the electromagnetic (EM) fields 
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and geometry are negligible. The simulation box is 75A x 40A large (corresponding to 30 01 1 501 grid points), terminated by 
UPML boundaries to suppress reflection from the boundaries f^. The device is illuminated from the bottom (Fig.[T]i, using a 
current source that generates a monochromatic plane wave that propagates in the vertical direction. 

In Fig.|2lleft) we show a comparison between the angular distribution of the transmitted intensity /(0;OOO) as obtained from 
the FDTD simulation (bullets) and Fraunhofer theory (solid line). Plotting 

_ /(0;OOO) -/(0;COO) -/(0;OCO) -/(0;OOC) +/(0;CCO) +/(0;COC) +/(0;OCC) 
^ 7(0 =0;OOO) ' 

as a function of (see Fig. |2jright)) clearly shows that the WH Eq. ([T]i of Refs. 111121, is in conflict with Maxwell's theory: 
E(0) takes values in the 0.5% range, much too large to be disposed of as numerical noise. Note that is obtained from data 
produced by seven different device configurations. 

Physically, the fact that ^ is related to the presence of wave amplitude in the vicinity of the surfaces of the scattering 
object (one, two, or three slit system), see for instance Fig. [TJright). These amplitudes are very sensitive to changes in the 
geometry of the device, in particular to the presence or absence of a sharp edge. Although these amplitudes themselves do 
not significantly contribute to the transmitted light in the forward direction, it is well-known that thek existence affects the 
transmission properties of the device as a whole oHi^. 



B. Wave decomposition 

The essence of a wave theory is that the whole system is described by one, and only one, wave function. Decomposing this 
wave function in various parts that are solutions of other problems and/or to attach physical relevance to parts of the wave is 
a potential source for incorrect conclusions and paradoxes. Even for one-and-the-same problem, the idea to think in terms of 
waves made up of other waves can lead to nonsensical conclusions, such as that part of a light pulse can travel at a superluminal 
velocity. Of course, we may express the wave field as a superposition of a complete set of basis functions, e.g. by Fourier 
decomposition, and this may be very useful to actually solve the mathematical problem (to a good approximation). However 
such decompositions are primarily convenient mathematical tricks which, in view of the fact that in principle any complete set 
of basis functions could be used, should not be over-interpreted as being physically relevant ITtjSfl. 

The WH Eq. ([T]i takes these ideas substantially further by decomposing the wave amplitude in three parts, each part describing 
the same system (a single slit) located at a different position in space. It is then conjectured that the wave amplitude for the whole 
system (three slits) is just the sum of these three different amplitudes. 

Advocates of the "physical" motivation for this conjecture might appeal to Feynman's path integral formulation of wave 
mechanics to justify their picture but in fact, one can see immediately from Feynman's path integral formalism that the WH 
Eq. ([T]i is not valid. 

We use the expression for the propagator of the electron as given by Feynman and assume that the particle proceeds from a 
location a and time f„ on one side of the screen with slits labeled 1,2,3 to a location b where a measurement is taken at time ti, 
on the other side. We assume that there exists some time tc between and f/, (as assumed by Feynman on p. 36, Ref. 9). The 
propagator for this process is denoted by Feynman as K{b,a) (@]. If we include for clarity the times then we would have to write 
K{{b,th), (fl,fa)). As pointed out by Feynman (Ref.|^ p. 57) we have a connection of this propagator to the wave function iff 
given by; 

\l/{b,tb)^K{{b,tb),ia,ta)). (5) 

Feynman represented the propagator Khy a path integral that sums over all possible space-time paths to go from ato b with the 
end-point times as given above. If we have an infinitely extended screen in between a,b with only slit 1 open, then all paths can 
only proceed through this one slit. We denote the wave function that is calculated for a path leading through a particular point 
xi of the slit at time fvi by y/f. Similarly for slits 2 and 3 open only we have ™d i/Zg respectively and the corresponding K's 
are calculated with Feynman paths that only go through slits 2 or 3 respectively. 

Had we chosen all three slits open, then Feynman's formalism insists that pathways going through multiple slits matter in 
general. Therefore, we would have to include paths through multiple slits in the path integral representation of K and we would 
obtain a corresponding I/J23- Thus, Feynman's quantum mechanics with all three slits open does contain an infinity of paths that 
go through multiple slits resulting in yf[jj^. However, none of the wave functions i/j or y/^ may contain any path through 
more than one slit because of the assumption that only one slit be open at a time. Therefore all the expressions involving these 
amplitudes do not contain multiple-slit path integrals and consequently do not contain all the paths that are required to compute 
^f[23■ '■^^ ^^^^ subsection, we illustrate the importance of "all" by solving the Maxwell equations for a minor variation of the 
three-sht experiment in which we block one slit. 
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FIG. 3: (Color online) Top and middle: Amplitudes of the Ex (left) and (right) components of the electric fields as obtained from a FDTD 
solution of Maxwell's equation for light incident on a metallic plate with two slits and a hole between the two slits. The incident wave is 
monochromatic and has wavelength A = 500nm. The slits are A wide, their centres being separated by 3A. The index of refraction of the 
41 -thick metallic plate (steel, colored black) is 2.29 + 2.61i. The holes are X wide and 2X deep. Bottom left: Normalized angular distribution 
of light transmitted by the devices extracted from the FDTD simulation data. Bullets: Simulation results for 1(6)/ 1(6 = 0) for the slit in the 
centre filled with material half-way from the bottom (see center panel); Crosses: Simulation results for I' (6) /I' (6 = 0) for the slit in the centre 
filled with material half-way from the top (see top panel); Dashed line: Guide to the eye. On the scale used, the two angular distributions 
cannot be distinguished. Right: Normalized difference between the angular distributions of the device with the hole in the bottom (top panel) 
and the top (middle panel): -/'(0))/max(/(e = 0),I'(6 = 0)). According to the WH of Refs. Q, 01, this difference should be zero. 



C. Three-slit device with blocked middle slit 

The geometry of the device that we consider is depicted in Fig. [3] together with the FDTD solution of the EM fields in the 
stationary state. We have taken the three-slit device used in Fig. |2] and blocked the middle slit by filling half of this slit with 
material (the same as used for other parts of the three-slit device), once from the top. Fig. [Stop), and once from the bottom. 
Fig. [3j middle). Comparing the FDTD solutions shown in Fig.jSjtop) and Fig. [3j middle), it is obvious to the eye that the wave 
amplitudes provide no support for the idea that these systems can be described by a wave going through one slit and another 
wave going through the other slit. The angular distributions for the two cases look very similar (see Fig. [3jbottom,left)) but 
differ on the one-percent level (see Fig.[3lbottom,right)). 
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III. THE WORKING HYPOTHESIS AS AN APPROXIMATION 



Having presented examples that clearly demonstrate that WH E g. dp does not hold in general, it is of interest to scrutinize the 
situations for which the WH Eq. ([T]l is a good approximation fTol 1 1 ih . As pointed out earlier, in general, interference between 
many different pathways is not simply the sum of the effects from all pairs of pathways. To establish nontrivial conditions 
under which it truly is a pairwise sum, we discard experiments for which the WH trivially holds, that is we discard experiments 
that exactly probe the interference of three waves, such as the extended Mach-Zehnder interferometer experiment described in 
Ref. J^and the class of statistical problems described by trichotomous variables considered in Ref.fTsl 

Let us (1) neglect the vector character of EM waves and (2) assume that the diffraction of the three-slit system is described 
by Fraunhofer diffraction theory. Then, for normal incidence, the angular distribution of light intensity produced by diffraction 
from shts is given by [4J 

IiB,s,d,N) = } . > : ■ a ] : (6) 

\ sm(;rdsm0) J \ Kss,md J 

where s and d are the dimensionless slit width and slit separation expressed in units of the wavelength A, respectively. Therefore, 
we have 

A(0) = /(0,i,c/,3)-2/(0,i,c/,2) -7(0, ,?,2£/,2) + 3/(0, i,<i,l) 

= [(1 +2cos2fl(i)2-8cos2fl(i-4cos2 2flc/ + 3] ^^^^^ =0, (7) 

where a = tt sin 0. Thus, in the Fraunhofer regime the WH Eq. ([T]) holds. 

It is not difficult to see that A(0) = is an accident rather than a general result by simply writing down the Maxwell curl 
equations L3,iJ 

e(r)^^M = VxH(r,f)-/(r,f) 
at 

M(r)^5|^ = VxE(r,0, (8) 

where the geometry of the device is accounted for by the permittivity e(r) and for simplicity, as is often done in optics 10], we 
may assume that the permeability (r) = 1 . 

Let us write e(r,000) for the permittivity of the three-slit geometry and E(r,f,000), H(r,f,000) for the corresponding 
solution of the Maxwell equations Eq. (|8]l. The WH Eq. ([T]l asserts that there should be a relation between (E(r,f,000), 
H(r,f ,000)) and (E(r,f ,COO), H(r,f ,COO)), (E(r,f,OCO), H(r,f ,OCO)), (E(r,f ,OCC), H(r,f,OCC)) but this assertion 
is absurd: There is no theorem in Maxwell's theory that relates the solutions for the case e(r,000) to solutions for the cases 
e(r,COO)... e(r,OCC). The Maxwell equations are linear equations with respect to the EM fields but solutions for different e's 
cannot simply be added. 

Of course, this general argument applies to the Schrodinger equation as well. For a particle moving in a potential, we have 

m^-^^ = fV + V(r)V(r,0- (9) 



m 



In essence, the WH Eq. ([T]) asserts that there is a relation between the solutions of four problems defined by the potential V (r) 
and three other potentials Vj{r) for j = 1,2,3. More specifically, it asserts that 

^(r,f) = ^i(r,f)+^2(r,0 + 'I'3(r,f), (10) 

and 

y(r)1'(r,f) = y(r)^i(r,f)+V(r)'P2(r,f)+^(r)^3(r,f) 

= yi(r)^i(r,f)+V2(r)1'2(r,0+V3(r)^3(r,f)- dD 

The authors could not think of a general physical situation that would result in Eq. (fTTT i. 

In summary, A(0) 7^ in experiments not exactly probing the interference of three waves, but A(0) = in experiments carried 
out in the Fraunhofer regime. In real laboratory experiments, such as the one reported in Ref. 0, it is very difficult, not to say 
impossible, to measure A(0) = 0. Therefore, in the next section we present a computer simulation study of this experiment 
resulting in a quantitative analysis of the applicability of the WH Eq. ([T]). 
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FIG. 4: (Color online) Two-dimensional representation of the experiment reported in Ref. 0]. The three slits at the top are 30^m wide, their 
centres being separated by 100/xm. The blocking mask at the bottom can have one, two or three slits, each slit being 60/im wide with its centre 
aligned with one of the slits in the top plate |2]. In the example shown, the middle slit of the blocking mask is closed (corresponding to the 
case OCO). The separation between the top plate and blocking mask is 50/xm. The index of refraction of the 25/im-thick material (colored 
black) is 2.29 + 2.61; (index of refraction of iron at 405nm). The wavelength of the incident light is 405nm. Also shown are the amplitudes 
of the Ex (left) and E, (right) components of the electric fields as obtained from a FDTD solution of Maxwell's equation for a monochromatic 
light source (not shown) illuminating the blocking mask. Note that the E^- and £;-components propagate in a very different manner. 



IV. COMPUTER SIMULATION OF THE EXPERIMENT REPORTED IN REF.Q 



The geometry of this device is depicted in Fig. |4] (see also Ref. 13), together with the stationary state FDTD solution of the 
Maxwell equations. In the simulation, the device is illuminated from the bottom (Fig.[T]), using a current source that generates a 
monochromatic plane wave that propagates in the vertical direction. The wavelength of the light, the dimension of the slits and 
their separation, blocking masks and material properties are taken from Ref. [2]. In view of the large (compared to wavelength) 
dimensions of the slits, to reduce the computational burden, we assume translational invariance in the direction along the long 
axis of the slits. This idealization of the real experiment does not affect the conclusions, on the contrary: It eliminates effects of 
the finite length of the slits. In all these simulations, the 81 mesh points per wavelength {X — 405 nm) were taken to ensure that 
the discretization errors of the EM fields and geometry are negligible. The simulation box of 820^m x 120/im (corresponding 
to 3 936 188 001 grid points) contains UPML layers to eliminate reflection from the boundaries [3]. Each calculation requires 
about 900GB of memory and took about 12 hours, using 8192 processors of the IBM BlueGene/P at the Jiilich Supercomputing 
Centre. 




9 (degrees) 9 (degrees) 



FIG. 5: (Color online) Left: Angular distribution of light transmitted by the system shown in Fig. |4] for the cases in which all slits are open 
(000: solid line), one sUt is closed (COO and OOC: dashed line, OCO: dash-dotted line), and two slits are closed (CCO,COC,OCC: double- 
dotted line), as obtained from FDTD simulations. Right: K" as a function of 9, as defined by Eq. l ll2t . According to the WH yj,|2j], this 
difference should be zero. 



Qualitatively, Fig. |4tleft) indicates that the x-component {Ex) of the EM-field propagates through the two layers of slits with 
very little diffraction from the top (= blocking) layer This is not the case for the z-component shown in Fig. SJright). In this 
idealized simulation setup, the amplitude of the y-component of the EM-field is zero. 

In Fig.|5];ieft) we present the results for the angular distribution of the seven cases (OOO, OOC, OCO, COO, OCC, COC, and 
CCO), extracted from seven FDTD simulations. From Fig.|3right), it is clear that k{Q), defined as [2] 

/(e;000) -/(e;COO) -/(e;OCO) -/(e;OOC) +/(e;CCO) +/(e;COC) +/(e;OCC) 

' " |/(0;OOC) -/(0;OCC) -/(0;COC)| + |/(0;COO) -/(0;COC) -/(0;CCO)| + |/(0;OCO) -/(0;OCC) -/(0;CCO) 

(12) 

is not identically zero, but of the order of 10^^. Note that Eq. (fT2b exactly corresponds to the expression for K defined in 
Ref. since in the idealization of the real experiment 7(0; CCC) = 0. In Ref. 2 it is reported that k-(0 = 0) = 0.0064 ±0.0120 
for measurements with single photons, k{Q = 0) = 0.0073 ±0.0018 for measurements with a laser source and a power meter 
for detection and k{Q = 0) = 0.0034 ± 0.0038 for measurements with a laser source attenuated to single-photon level and a 
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silicon avalanche photodiode for detection. The upper bound for K at several detector positions given by the experiment 1*5] is 
k{9) < 10"^. We find for the idealized version of the experiment k{9 = 0) = 4 x 10"^ and k{9) <1 x 10^^, a factor 100 smaller 
than the values measured in the experiment. Note that the experimental and simulated values for K are very small because the 
experiment is carried out in a regime in which scalar Fraunhofer theory works well, as can be expected from the dimensions of 
the slits and slit separations of the device. 

V. DISCUSSION 

A necessary condition for a mathematical model to give a logically consistent description of the experimental facts is that there 
is one-to-one correspondence between the symbols in the mathematical description and the actual experimental configurations. 
When applied to the three-slit experiment in which one slit or two slits may be closed, the argument that leads from Eq. ([T]) to 
Eq. (|2|i is false because there is no such correspondence. 

If in Eq. ^ is to represent the amplitude of the wave emanating from the yth slit with all other slits closed, the WH should 
be written as 

/(r,000) = Iv/^(r,OCC) + v/^(r,COC) + v/^(r,CCO)p, (13) 

that is, we should label the y/s such that there can be no doubt about the experiment that they describe. This notation establishes 
the necessary one-to-one correspondence between the mathematical description (the \//'s) of the particular experiment (labeled 
by OCC, etc.). Now, we have 

/(r,000) = |v/(r,OCC) + v/(r,COC)p+|V/(r,OCC) + v/(r,CCO)|2 
+ |V/(r,COC) + v/(r,CCO)|2 

-|V/-(r,OCC)p-|V/(r,CCO)p- |V/(r,CCO)p. (14) 

At this point, it is simply impossible to bring Eq. (fl4] i into the form Eq. (|2]i without making the assumption that 

V/(r,OOC) = i/(r,OCC) + v/(r,COC), 
V/(r,OCO) = i/(r,OCC) + v/(r,CCO), 

V/(r,COO) = i/(r,COC) + v/(r,CCO). (15) 

If we accept this assumption, we recover Eq. (O. However, the assumption expressed by Eq. ( fTsT i cannot be justified from 
general principles of quantum theory or Maxwell's theory: The only way to "justify" Eq. (fTSl l is to "forget" that the i/'s are 
labeled by the type of experiment (e.g. OCC) they describe. For a discussion of this point in the case of a two-slit experiment, 
seeRefs.lllll. 

In other words, accepting Eq. ( fTSl l destroys the one-to-one correspondence between the symbols in the mathematical theory 
and the different experimental configurations, opening the route to conclusions that cannot be derived from the theory proper. 
Hence, if A(r) for a three-slit experiment, one cannot conclude that Born's rule does not strictly hold. 

From the above it is clear that in general (excluding special cases such as the Fraunhofer regime), the three-slit interference 
pattern cannot be obtained using a combination of single-slit and two-slit devices. In order to measure the three-slit interference 
pattern, a three-slit device with all three slits open is required. Similarly, for the measurement of the two-slit interference pattern 
a two-slit device is required; a combination measurements using two single-slit devices is not sufficient. A similar issue was 
raised in Ref. 16 concerning the measurement of arbitrary Hermitian operators acting on the Hilbert space of a spin-5 particle 
using suitable generalized Stern-Gerlach apparatuses. The standard Stern-Gerlach apparatus with an inhomogeneous magnetic 
field whose direction is constant, but whose magnitude depends on the position, can measure any Hermitian operator of a spin- 
1/2 system. This is however not the case for spin-5 systems with 5 > 1/2 (l^. In order to measure all spin-5 operators a 
generahzed Stern-Gerlach apparatus, using both electric and magnetic fields, is required l,16J . 

VI. SUMMARY 

The results of this paper can be summarized as follows: 

1. The three-slit interference, as obtained from explicit solutions of Maxwell's equations for realistic models of three-slit 
devices, is nonzero. 

2. The hypothesis flSlj that the three-slit interference Eq. (|3]l is zero is false because it requires dropping the one-to-one 
correspondence between the symbols in the mathematical theory and the different experimental configurations opening a 
route to conclusions that cannot be derived from the theory proper. 
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3. Although not holding in general, the hypothesis that the three-slit interference Eq. Q should be zero is a good approx- 
imation in experiments carried out in the Fraunhofer regime. The experiment reported in Ref. 2i is carried out in this 
regime and bounds the magnitude of the three-slit interference term to less than 10^^ of the expected two-slit interference 
at several detector positions |12[]. By explicit solution of the Maxwell equations for an idealized version of the three-slit 
experiment used in Ref. |2| we provide a quantitative analysis of the approximative character of the hypothesis that the 
three-slit interference Eq. (O is zero. We find that the magnitude of the three-slit interference term is several orders of 
magnitude smaller than the upper bound found in the experiment 101 • 
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